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We study an observational method to analyze non-Gaussianity of a gravitational wave (GW) 
background made by superposition of weak burst signals. The proposed method is based on fourth- 
order correlations of data from four detectors, and might be useful to discriminate the origin of 
a GW background. With a formulation newly developed to discuss geometrical aspects of the 
correlations, it is found that the method provides us with linear combinations of two interesting 
parameters, I2 and V2 defined by the Stokes parameters of individual GW burst signals. We also 
evaluate sensitivities of specific detector networks to these parameters. 
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^ ■ I. INTRODUCTION 

In the last decade, sensitivities of gravitational wave (GW) detectors have been dramatically improved [H,[1,[I|. It is 
expected that a similar trend would continue in the next decade, and we will soon detect GWs directly. Then a totally 
new branch of astronomy will be opened with various observational targets from astrophysics to fundamental physics. 
One of the most interesting targets of GW observation is a stochastic GW background [1, 0, H, @|- For example, we 
■ might obtain crucial information of physics at an extremely high-energy scale by analyzing a GW background generated 
in the early universe. Around / ^lOOHz, the measurement sensitivity of the Laser Interferometer Gravitational Wave 
Observatory (LIGO) to a GW background is about to surpass the level OgwC/) ~ 10~ 5 that can be indirectly 
5^ ' constrained by the observed abundances of light elements composed at the epoch of nucleosynthesis [1] (see also 0]). 
Here the function flcw(f) is the energy density of a GW background per logarithmic frequency interval normalized by 
the critical density of the universe [J] . In order to efficiently extract information from a GW background, it is essential 
I , to quantify it with various measures in addition to the traditional spectrum ^lGw(f) (see e.g. Q for anisotropies and 
£> ' @, 03 for asymmetry of right- and left-handed GWs). 

In a recent paper [ll| (hereafter Paper 1), the author discussed detectability of non-Gaussianity of a GW background 
caused by short-duration bursts whose signals are individually weak and undetectable (see also [HI, [TH, 0, EH ) • The 
origin of the non-Gaussianity here is the discreteness and finiteness of the numbers of the burst events. The proposed 
method might be useful to discriminate the genesis of a background e.g. whether it is smooth enough to be consistent 
with that generated during an inflation epoch. In addition, the author pointed out that by a statistical amplification 
of weak bursts signals, we might estimate their basic characters such as rate and amplitude of bursts and their typical 
duration. In this follow-up paper, we extend the previous study with fixing some overlooked points in Paper 1, 
and also provide a new formulation for quantitatively discussing the non-Gaussianity measurement in relation to the 
geometry of a detector network and polarization properties of burst GW emissions. 

This paper is organized as follows; in section II, using a simplified model, we explain our basic approach for 
the non-Gaussianity measurement based on fourth-order correlations, and discuss application of this approach to 
observational studies on GW backgrounds. In section III, we analyze how the non-Gaussianity measurement depends 
on the geometry of a detector network. We define the generalized overlap reduction functions £77 and £yy that 
characterize sensitivity of the measurement for a given four-detector network. Then, in section IV, the generalized 
overlap reduction functions are evaluated for two specific networks of GW detectors both on the Earth and in space. In 
section V, we calculate expressions for estimating signal-to-noise ratios of our non-Gaussianity measurement. Section 
VI is a brief summary of this paper. Two appendixes are added. In appendix A, we examine correlations of Fourier 
modes transformed with a short time duration. This appendix is technical, but would be useful to understand some 
of basic properties in this paper. In appendix B, we comment on the overlooked points in Paper 1. 



II. DETECTING NON-GAUSSIANITY 



A. underling approach 



In this paper, we mainly study a GW background made by a superposition of burst-like signals, and discuss how 
to analyze its non-Gaussianity. As a preliminary set-up, we first outline the underlying approach using a simplified 
model. We consider four detectors a,b,c, and d with representing their data u a i,uu,u C i and Udi respectively. Here 
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i is the suffix for the data sequence (e.g. Fourier modes). But, since the important issue in this subsection is the 
inter-detector correlation structure within a single i and not the correlation between different i, we omit the suffix i 
below. 

We assume that data uj (J — a, b, c, d) are superpositions of the following three ingredients; (i) detector noises 
nj with no correlation between other detectors, (ii) a Gaussian signal g (mimicking a smooth inflation-type GW 
background) and (iii) a group of independent burst signals Vj . Here I is the number of bursts in the data and a 
random variable. Then the four data are written as 

i i i i 

u a = n a +g+}^Vj, u b = n b + g + ^Vj, u c = n c + g +}^vj, u d = n d + g +}^vj. (f) 

3 3 3 3 

The latter two ingredients g and . Vj are set to be identical for all the detectors, and we dropped the label J for 
these two components. This is just for simplicity and is not essential for our demonstration of extracting non-Gaussian 
signature. In the next subsection, we deal with the differences between detectors taking into account their responses 
to incident burst signals. 

We assume that the number of bursts / obeys the Poisson statistics whose probability distribution function P(l) is 
solely determined by the averaged event number q as 

P(l) = (2) 
Using this functional form, we can derive the following identities for moments 

oo oo oo 

5Xo = i, 5Xoi = «, 5xo*o -i)=? 2 - (3) 

1=0 1=0 1=0 

The first expression shows the proper normalization of the probability distribution function, while the second one is 
the definition of the averaged event rate q. 

So far we have discussed the data uj in somewhat abstract manner. In the standard correlation analysis of a 
GW background, the data uj are the Fourier modes transformed from data streams acquired in the time domain. 
Therefore, we treat our data uj as complex numbers, and define moments required for evaluating correlations of the 
data uj. We represent moments of the Gaussian signal g as 

G 1 = {gg*), G 2 = {gg), (4) 

where the notation (■ ■ •) represents to take an ensemble average. Even if the original (e.g. time domain) data stream 
is a real and Gaussian variable, the amplitudes of the real and imaginary parts of its Fourier mode g can have different 
expectation values for a transformation with a short time segment (see Appendix A). This is the reason we introduced 
two moments G\ and G 2 . In the same manner, we denote the moments of individual signal Vj as 

D 1 = (v j v*), D 2 = (v jVj ), D 3 = (( Vj v*f). (5) 

Note that the above correlations can be estimated by taking the averages with respect to the dropped suffix i for 
Fourier modes, as explicitly discussed in the next subsection (e.g. eqs. (IT51) and (ffU|) ). 
We can now write down the correlation of two data a 7^ b as 

(u a ul) = qDi + d, (u a u b ) = qD 2 + G 2 - (6) 

Here, we used the assumption that detector noises n.j are statistically independent, and separate burst signals Vi and 
Vj (i j) are uncorrelated. In the same manner as the second-order moments, the forth-order moment (u a UbU*u^) 
is evaluated as 

(u a u b u* c u* d ) = qD 3 + q(4D 1 G 1 + D 2 G* 2 + D* 2 G 2 ) + q 2 (2D\ + D 2 D 2 ) + 2G\ + G 2 G* 2 , (7) 
and, using eqs.([6]) and ||7J), we obtain 

JC = (u a u b u* c u* d ) - (u a u*) (u b u* d ) - (u a u b ) (u*u d ) - (u a u* d ) (u b u* c ) = q{(vjv*) 2 ) = qD 3 . (8) 



Note that this combination is proportional to the averaged event number q and does not have contribution from the 
cross terms. We can extract a quantity that is purely arisen by discreetness of the underlying signals and closely 
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related to their non-Gaussianity. Actually, the combination K. is essentially the same as the Kurtosis parameter, a 
well known measure in astrophysics to characterize non-Gaussianity (see e.g. [l6| for its application to the large-scale 
structure in the universe). For Gaussian variables Xi (i = 1,2,3 and 4), we have a simple identity (X1X2X3X4) = 
(X1X2) (X3X4) + (X1X3) (X2X4) + (X1X4) (X2X3), and we can easily confirm that the parameter K, vanishes for 
Gaussian variables. Indeed the right-hand-side of eq.© does not depend on G\ and G2. A third-order moment called 
skewness is often used to characterize asymmetry of a probability distribution function for a scalar-type quantity 
around its mean. However, since GWs are tensor quantities with no preferred signs, it is irrelevant to analyze the 
skewness parameter here for characterizing non-Gaussianity of a GW background. 

We should also comment on the relation between the parameter K, and the central limit theorem. Let us 
consider a situation where we increase the number of burst q, but fix the total power of the bursts • Vj by 

(y~]j Vj)(^2j v j)*\ = qDi = const and also keep the ratio D$/ D\ = const for the individual burst. The parameter 

K = qD?, = (D3/ Df)(qDi) 2 q~ l oc q^ 1 asymptotically approaches to zero for q — > 00, and the data (u a ,Ub,u c , Ud) 
become more Gaussian-like. This matches with our expectation from the central limit theorem. 



B. GW observation with multiple detectors 



In the previous subsection, we provide a simple demonstration for extracting non-Gaussianity signature caused 
by discreteness of burst signals. Here, we specifically discuss application of the approach to a GW background. As 
there are some gaps between the settings in the previous subsection and those in other parts of this paper, we do not 
straightforwardly use the notations of variables defined in the previous subsection, but we introduce new (and mostly 
distinct) ones that will be used hereafter. 

We assume to have multiple detectors (labels J = a, 6, • • • ) that commonly have optimal sensitivity around a 
frequency / with a bandwidth A/ ~ /, as expected for typical laser interferometers. We consider a signal analysis 
in the optimal band, and neglect details of frequency dependence (e.g. replacing an integral J (• ■ • )df with a product 
(•••)/ x Af ~ (•••)/ x /)■ I n practice, this situation is approximately realized by applying a band-pass filter. 

We model the data stream sj(t) of a detector J in terms of a GW signal Hj and a detector noise nj as 

Sj (t)=Hj(t)+nj(t). (9) 

For analyzing the latter nj, it is advantageous to work in the Fourier space. We decompose the data streams (total 
duration T b s ) into short segments of a given duration T Beg (> f~p t ), and attach a label M(= 1, • • • ,T bs/T seg ) for 
each segment. Then we take Fourier transformations at discretized frequencies / = NT~ eg (N: integer) as follows 

sjM(.f) = / e 2mft sj(t)dt = H JM (f) + n JM (f), (10) 

J(M-l)T 3eg 

with 

HjmU) = / e^Hj^dt, njutf) = / e^ iSt nj(t)dt. (11) 

J{M-l)T teg J(M-l)T SBg 

The number of relevant Fourier modes in a segment is ~ T seg Af ~ T seg f. 

We assume that the detector noises njM are stationary, Gaussian, independent, and have identical spectrum Sjv(/)- 
Then the covariance matrix for the detector noises is given by 1 

(nj M (f)n K L(fT) ^ \5 ML 5 JK 5 ff ,T seg S N {f). (12) 

Among the statistical assumptions about the detector noises, independence between detectors is the critical one for 
structure of our approach, but other ones (including Gaussianity) would only modify statistical significance of data 
analysis [ill ]. 

The total number of Fourier modes in observational time T b s is given by 

Nt ~ ^ x T seg Af = T obs Af. (13) 

-l seg 



1 In reality, especially for short segment length T seg Af ~ 1, there is a weak correlation between nearby modes. 
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Since we will take statistical averages of these Fourier modes with neglecting frequency dependence, we use a running 
index i(= 1, • • • , N t ) for the modes without applying the double decomposition by M and / as in eq. pH)) . One of the 
reasons behind this prescription is that, as we see later, it is advantageous to take a short segment time T seg ~ (A/) -1 
for increasing the signal-to-noise ratio of the non-Gaussianity measurement. 

In the previous subsection, we find that the smooth Gaussian component g in eq.(l) turned out to be irrelevant for 
the Kurtosis parameter. Therefore, in the followings, we put aside such a contribution, and only consider the burst 
signals for GW signals 

U 

Hji = Y,Rjii, (14) 

3 

where U is the total number of bursts in the time segment relevant for the Fourier mode i, and R.jij is response of a 
detector J to a burst characterized by the labels i (for Fourier modes) and j (for bursts in a mode i). 
We directly apply results in the previous subsection for our data now given by 

s a i = H ai + n ai , Sbi = Hbi + ribi, s C i = H ci + n ci , Sdi = H di + n di . (15) 

We can basically identify the current data sji in eq.QISp with the previous one uj in eq.(l) (J: the label for detectors). 
But, unlike the burst signals Uj in eq.(l), we keep the label J for the burst component Rjij to include geometrical 
dependence of a detector network that will be studied in the next section. From the data Sji we take the averages 
over the Fourier modes 

1 N t N t 
C 2 \ab = Jjr 2J Saisti, C 2 2ab = JT S ai s bi- (16) 

These contain statistical fluctuations around their expectation values. In the present setting, we have the expectation 
values 

(C 2 lafc) = (SaiS* bi ) = (HaiHfo) = q (R a Rt) , (C 2 2ab) = (SaiSbi) = (H al H bl ) = q (R a Rb) (17) 

with the average number of events q = rT seg [r: event rate) in a segment T seg . Here we used statistical independence 
of detector noises, neglected frequency dependence, and also omitted the label j for bursts (e.g. replacing (^RaijRlij 
with (R a Rl)). We use the summations C'2iab and C22ab as the estimators for the expectation values (sats^) and 

(SaiSbi) as follows 

G 2 \ab — > (SaiSbi) = q (R a Rb) , C*22afc — > (s a iSbi) = 9 (R a Rb) ■ (18) 

The fluctuations of the summations C2i ab and C 2 2ab around their expectation values would be evaluated in subsection 
V.B. Similarly, we define the summation 



1 Nt 

Ciabcd = TT ^2(SaiSbiS* ci S* di ) (19) 



N t 

i 

for the estimator of the forth-order moment ((saiSbiS^s*^). We then introduce the combination 

K = Ciabcd — C2\acC 2 \bd ~ C2\adC2\bc ~ C22abC 2 2cd (20) 

as an estimator of the kurtosis corresponds to eq. (JSJ) 

(K) -> (s ai s bi s* ci s* di ) - (s ai s* cl ) (s bl s* di ) - (saiSu) (s* ci s* di ) - (s ai s di ) (s H s* ci ) = q (R a R b R*R d ) ■ (21) 

Therefore our observational target K is characterized by a correlation of four detectors (R a RbR* c R*d) defined by an 
ensemble average of single burst event. 



III. FORMULATION FOR GENERALIZED REDUCTION FUNCTIONS 



In this section, we discuss dependence of the expectation value (R a RbR* c R* d ) on the geometry of a detector network 
and the polarization state of incident burst GWs. Since the extra-Galactic burst sources would have random directions 
and orientations for realistic astrophysical models, we need to deal with many averaging operations with respect to 
angular parameters of the detector-sources system. 
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A. beam pattern functions of detectors 

We assume that the individual burst duration Td is much shorter than the time scale for transformation of the 
detector network {e.g. ~1 day for ground based detectors, ~1 year for typical space interferometers). Thus we fix the 
network configuration and introduce a fixed spherical coordinate system that is attached to the detector network (D- 
system, see figure 1). We denote the source direction of the GW signal by h (equivalently, the propagation direction: 
— n) whose explicit form is given by 2 

n = (cos (j>D sin^D, sin^u sin#D, cos 9p). (22) 

Since we introduce many angular parameters for the geometry of the network and burst GW sources, we also use the 
simplified notation "D" to represent the two angular parameters (9d,4>d) related to the detector network. In order 
to discuss GWs from the direction n, we define two unit vectors eg D and e^ D that are normal to the source direction 
n and given by 

eg D = (cos 6d cos <pr>, cos 9d sin 4>n, — sin^o), e l f >D = (— sa\.(pD, cos <pr>, 0). (23) 

Then the bases for transverse-traceless tensor e p (P = +, x) associated with the propagation direction n are given 
as (see e.g. 0, @|) 

e+ = eg D <g> eg D - e 4>D ® e <j)D , e x = eg D ® e^ D + e cf>D ® eg D . (24) 

In this paper we only deal with simple L-shaped interferometers. Response of such a detector is characterized by a 
tensor dj given by directions iij and vj of its two arms (see also [I3])- This tensor is traceless and written by 

dj = (uj (Siij-vj® Vj)/2. (25) 

The beam pattern functions Fj represent sensitivities of the detector to the two polarization modes defined in eq. (f24"|) . 
They are formally written with the tensor dj as 

Ff(9 D , fo) = dj : e p (n) = ^ d. Hj e p , (26) 

ij 

with two polarization states P = +, x. Note that Fj~ has even parity, while Fj has odd parity. These properties will 
become useful for geometrical interpretation of our results derived later. 

B. response of detector to GWs 

Next we characterize GW emission from a burst source. In addition to the D coordinate defined in previous 
subsection, we introduce another coordinate system S = (9s, 4>s) attached to the source (see figure 1). We also define 
two polarization bases (plus (p) eP and and cross (c) e c modes) 3 associated with the source coordinate similar to 
eq. ([24l . With these bases we decompose the polarization patterns of the burst GWs and represent their Fourier modes 
by (h p (f), h c (f)) that are defined at the origin of the detector network coordinate D. 

Now we discuss the correspondence of two polarization bases (e + ,e x ) and (e p ,e c ). The polarization angle ip 
characterizes the rotation angle between two coordinate systems D and S around the direction n (see figure 1 and 
18]), and the decomposition of the burst signals in the D system (h + ,h x ) is given by 

h+ = (hp cos 2ip + h c sin 2tp), h x = (— hp sin 2ip + h c cos 2ip). (27) 

Then the response of the detector J to the incoming burst GW is given by 

Rj = (F+h+ + F?h x ) exp[-2mfxj ■ h] , (28) 



2 We attach a hat " for an unit vector. 

3 Here we use the labels p and c for plus and cross modes to distinguish the + and X modes defined in the detector network coordinate. 
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where we formally added the phase factor exp[— 2mfxj ■ n] induced by the position of the detector relative to the 
origin of the D coordinate. This prescription is valid, when the light-travel time between detectors are much smaller 
than the segment time T seg for Fourier transformation. With eqs. ([27|) and (|28|) . the response Rj is written as 

Rj = T v jh p + T c jh c . (29) 

The functions J-j' c depend on the three angular parameters D = (9d, 4>d) and tp. The information of the phase factor 
exp[— 2wifxj ■ h] is included in them as follows 

T p j = (F+cos(2^))+Fj < sin(2V'))exp[-27ri/Kj-n], (30) 
JTj = (~Fj sin(2^) + Fj cos(2^>)) exp[-2wifxj • n]. (31) 

As mentioned before, it is reasonable to assume that orientations and directions of the extra-Galactic burst sources 
are randomly distributed. In other words, the source distribution is isotropic with no global handedness. In this 
situation, we define the following three averaging operators for the direction angles D = (6ut4>d), the orientation 
angles S = (6 s, <j>s) and the polarization angle ip; 4 

[Z(D,S,tP)] D = — dcb D d6 D sme D Z(D,S,iP), (32) 
47r Jo Jo 



[Z(D,S,i>)] s = —j d4>sj o d6 s sm6 s Z(D,S,iP), 



(33) 



[Z{D, S, V)U = ~ I WZ(D, S, V). (34) 



We also use the simplified notations such as [Z]ds = [[^]zj]s = [[-^]s]r>- 



C. overlap reduction functions 



Here we calculate the geometrical averages [RaR^DStp and [R a RbR* c R*d\DSi> to evaluate the expectation values 
(R a Rl) and (R a R\yR* c R*i) in eqs. (|T7l) and (f2~Tj) . For standard correlation analysis of a GW background, we take the 
combination R a R^ using two detectors a and b. Its angular average is formally evaluated as (9l. [To| 

[R a Rl]DSi> = l(llab[l]s + lVab[V] S ), (35) 

where the Stokes parameters / and V are defined by [lj| 

I=\h p \ 2 + \h c \ 2 , V = i(h p h*-h c h* p ). (36) 

The parameter J(> 0) represents the total intensity of two polarization modes and has even parity. The Stokes V 
parameter is related to the circular polarization of waves and has odd parity. The V parameter can be both positive 
and negative depending on the relative amplitude of right-handed waves (Hr = (h p + ih c )/y/2) and left-handed waves 
(fiL = (hp — ih c )/^/2), and we have a relation \V\ < I with equality only for 100%-circularly polarized waves. Note 
that these two real parameters I = |/i_r| 2 + |/il| 2 and V = — |/ia| 2 + |/il| 2 are invariant (spin-0) with respect to rotation 
of two polarization bases around the propagation direction, while the combinations Iir and themselves change as 
spin±2 quantities. Considering the fact that we take the averages for the angle ip corresponding to the rotation, it is 
reasonable that our expression should be a simple combination of the spin-0 parameters [I]s and [V]s- 

In this paper we consider GW sources with no preferred handedness, and we have the identity \V]s = 0. In other 
words, even if individual sources have circular polarization modes V (e.g. inspiral binaries), its ensemble average 
should vanish. In eq. (|35[) . the two functions 7/ and jv are overlap reduction functions defined by [5l Rj fTol] 

Ha b = \[TZ3V+?ZJT\d+ = +F Q x F b X)e- 2 ^-^)-"/ c ] D , (37) 



4 The angle i/> is introduced to specify the orientation of the source frame relative to the detector frame. Since GWs are spin-2 quantities, 
we can limit the integral of if) in the range [0, 7r], due to the apparent identity between -0 and ip + it. 
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Source frame 




FIG. 1: In the detector network (D) frame, the direction of a GW source n is characterized by two angles D — (0d,4>d)- In 
the source (S) frame, the direction of the detector network is given by S = (8s, 4*s)- On the plane normal to the propagation 
direction — n, projections of two z-axes (8d = and 6s = 0) are off-set by the polarization angle ip. These five angular 
parameters fix the relative configuration of two frames other than the distance between them. 



lv ab = KH*)\m> = ~ F-F+)e-^^- x ^l c ] D . (38) 

These characterize correlated responses of two detectors a and b to incoming GWs for averaged source configurations. 
To simplify our notations, we hereafter omit the subscript ab for the overlap reduction functions. With eq.([35p and 
the identity [V]s = for our source model, the angular average of the product R a Rl is given by 

[R a Rb]DSi> — -=~h> (39) 
5 

where we defined 

h = [I]s. (40) 

So far, we have studied the averaging operations with respect to the geometrical parameters of sources relative to 
a fixed detector network. In addition to these geometrical ones, we need to deal with averaging for the intrinsic wave 
pattern of bursts, including their distance distribution. But they are beyond scope of this paper. We rather assume 
that the parameter I\ implicitly encompasses these additional averaging, and simply put 

(R a R* b ) = ^-h. (41) 

Now we extend our analysis to the four-point combination (R a RbR*R*i) ■ After simple calculation, we find that the 
geometrical average [R a RbR*R*^\DS4i becomes a linear combination of the three parameters [I 2 ]s, [V 2 ]s and [IV] s as 
follows 



[R a R b R* c R* d ]nsi> = |r {Oi[I 2 ]s + Cvv[V 2 ] s + (iv[IV]s) 



(42) 
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Note that only the three parameters I 2 , V 2 and IV are the spin-0 quantities made from the relevant fourth-order 
moments of h p and h c . In eq. (|42p we put the pre-factor 2/25 in order to simplify the expression given later in 
subsection V.B. The coefficients Qi, (jv and (vv are generalization of the overlap reduction functions for the four 
point correlations and defined by 

95 

Cii = ^l(3F*F*F*F*+F+F+F*F*+F+F b *F+F*+F:F+F+F*+F+F*F*F d + 

+F*F+F?F+ + F*F*F+F+ + 3F+F+F+F+) exp[-27ri/(s a +x b -x c - x d ) ■ n}] D , (43) 



95 

Cvv = --[(F?F b *F c *F* + 3F+F+F?F* F+F*F+F} - F*F+F+F* - F+F*F*F+ 

-F*F+F*F+ + 3F*F*F+F+ + F+F+F+F+) exp[-27ri/(s +x b -x c - x d ) ■ h]] D , (44) 

95 

Civ = y l(F+F*F*F* + F?F+F?F* F*F*F+F* + F+F+F+F* F* F* F c * F+ 

+F+F+F C *F+ - F+F*F+F+ - F*F+F+F+) cxp[-2mf(x a + x b - x c - x d ) ■ n]) D . (45) 

In this paper, we call them as the generalized overlap reduction functions. Since we do not assume a global handedness, 
the average [IV]s should vanish due to a parity reason. But, here, it is important to note that our result in eq. (j4"2"|) 
depends on the quantity [V 2 ]s > that should not vanish for a burst model with V ^ before the ensemble average. 
In contrast to [IV] s or [V]s, the cancellation between right- and left-handed modes does not occur for the parameter 
[y 2 ]g. This, in principle, allows us to statistically study the circular polarization state of the burst by studying 
the non-Gaussianity of their background. As an example, we examined the ratio [T^ 2 ]s/[/ 2 ]s for a quadrupole-type 
emission pattern (h c ,h p ) cx [(1 + cos 2 0s)/2, ±i cos#s], and obtained the result [^ 2 ]s/[^ 2 ]s = 69/71. 
Similar to the definition of the parameter I\ = [I]s, we define the two parameters 1% and V% by 

h = [I 2 ]s, V 2 = [V 2 ] s , (46) 

including implicit averaging operations other than the geometrical ones. Then the fourth correlation (R a Ri,R*R2) is 
written as 

(R a R b R* c R* d ) = ^(hCii + V 2 (vv)- (47) 
To simplify some of expressions derived later, we also introduce a non-dimensional quantity W of order unity by 



W = i (hCii + V 2 Cvv) , (48) 

J 1 



and then we have 



(R a R b R* c R d ) = ^I 2 W. (49) 



IV. EXAMPLES OF THE OVERLAP REDUCTION FUNCTIONS Crz AND (vv 

Our fiducial target in this paper is a GW background made by a superposition of short GW bursts. As we see 
in the next section, it is preferable to take a short segment time T seg when analyzing such a background with a 
Fourier transformation. In this case, we cannot naively apply the phase factor such as exp[— 2irif(x a — x b ) ■ h/c] to 
evaluate the overlap reduction functions. This is because the coherent phase structure is not clear-cut for a Fourier 
transformation with a short time segment. Therefore we simplify our analysis below by working in the long- wave 
limit (namely f\(x a — x b )\/c ^ 1), and neglect the phase factors. Under this limit, the parity structure of the overlap 
reduction functions 7/, jv, (11 Cvv and (iv in eqs. ([37| (|38p ([43|) and (|44|) become particularly simple from the basic 
parity correspondences F^ — > Fj' and Fj — > ~~Fj f° r the beam pattern functions. As a result, we can easily show the 
identities ~fv — Civ = for the odd parity modes. In addition, the generalized overlap reduction function becomes 
invariant with respect to replacement of the detector indexes, due to its apparent symmetry in eq.(|43p. For example, 
we have [R a R b R* c R d ]DStp = [RaR b RcRd\DS'i>- I n this section, with the long-wave approximation, we evaluate the 
generalized overlap reduction functions Qi and (,vv for specific networks of interferometers both in the space and on 
the Earth. 



9 



A. BBO/DECIGO type detectors 

Our first example is the Big Bang Observer (BBO) [2(| and the Deci-hertz Interferometer Gravitational Wave 
Observatory (DECIGO) [2l|. They are future plans of GW observation in space with the optimal band around 
0.1-lHz. One of their main goals is to directly study a GW background from early universe by correlation analysis. 
However, it was pointed out that burst GWs from supernovae of population III stars might be a strong foreground 
for detecting a GW background from the early universe [13, [HI]. Therefore, the method proposed in this paper 
could be particularly useful for BBO/DECIGO, and might enable us to discriminate whether a detected background 
is consistent with a smooth (i.e. composed by many degree of freedom) Gaussian-like one. In this subsection, we 
briefly discuss the detector configurations of the proposed missions and evaluate their generalized overlap reduction 
functions. We will return to the astronomical aspects of the supernova background in subsection V.C. 

With BBO (and similarly for DECIGO), two sets of equilateral-triangle-shaped system (like LISA [24]) Y\ and Y~2 
would be configured to realize a good sensitivity to a background with a large overlap 77 ~ 1. We can make two 
effective L-shaped detectors (A\,E{) and (A2, E2) from the units Y\ and I2 respectively (see figure 2). For these four 
effective detectors, analysis under the low-frequency-approximation will be fairly well around their optimal frequencies 
/ = 0.1-lHz. 

These four detectors are placed on the same plane (see figure 2), and responses to GWs can be regarded as 
identical for A\ — A2 and E± — E2 pairs ("/ia 1 a 2 — liE t E 2 — 1) m the low frequency approximation. Meanwhile, 
in geometrical point of view, the two-detector pair Ai — E\ (same for A2 — E2) is misaligned by 45° on the plane 
[liA 1 E 1 = liA 2 E 2 =0)- As for the detector noises of the four data, we assume that their spectra are the same 
and their correlation is negligible. This is because of the independence of the noises of two systems Yi and Y2 and 
the geometrical symmetry within each triangle [25j. 5 These noise properties are highly preferable to measure the 
non-Gaussian parameter K with the four detectors A\, A2, E\ and E2. 

Due to the geometrical simplicity of the network, we can analytically calculate the functions (u and CvVi an d 
obtain 

Cn = I- (50) 

This result does not depend on the positions of two conjugates * among the four data, as commented earlier. On the 
other hand, we have the overlap functions for circular polarization modes as 6 

Cvv = {A 1 A 2 E* X E*2) (51) 

and 

QK 

Cvv = [A 1 AIE 1 E* 2 ). (52) 

Using the difference of the function Cyv as in eqs.(|5ip and l|52p . we can, in principle, measure the two parameters I2 
and V2 separately from their linear combinations. The basic prescription for such separation was proposed in [9j to 
deal with the circular polarization mode in standard (two-point) correlation analysis (see also [261]). 



B. Four detector network on the Earth 



We next evaluate the generalized overlap reduction functions Qi and (vv for a detector network composed by the 
on-going and planned ground-based interferometers listed in Table 1. For detectors widely separated on the Earth, the 
characteristic frequency below which we can apply the low-frequency-approximation is given by ~ c/(2ttRe) ~ 10Hz 
(Re- radius of the Earth ~ 6000km). This frequency is much smaller than the optimal frequency ~ 100Hz of the initial 



5 Roughly speaking, (Ai,Ei) are linear combinations of three basic data (s&<,j/i> «i) obtained at three vertexes using adjacent arms. 
They are written as Ai = Xl "^ % and Ei = Xl+v J^ z ' . The noises (n x i,n y i,n z i) of the basic three data apparently have correlation. 

But, if they are completely symmetric (namely, (n x iTi*^ = (jhii' n yi) = ( nzin ti) ana - {^ nxin yiJ = ( n y in ti) = ( nzin xi)^> we nave 
(nAin*Ei) = due to symmetric cancellations [25| |. In actual observation we need to carefully study the potential effects caused by 
residual noise correlations. 

6 Here, we can switch A\ <-> A2 and E\ <-+ E2 for getting the same results. 
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A 2 




FIG. 2: Schematic picture for the four effective detectors Ai, Ei, A2 and E2 made from two triangle units of BBO. In the 
long-wave approximation, the positional differences between detectors can be neglected and responses to GWs can be identified 
for the Ai — ^2-pair and Ei — i^-pair. Two detectors A1E1 has an off-set angle 45° (same for A2E2), and their noises are 
uncorrelated. 

detector 8 (f> a 

LCGT (C) 53.6 137.3 70.0 

LIGO Hanford (H) 43.5 -119.4 171.8 
LIGO Livingston (L) 59.4 -90.8 243.0 
Virgo (V) 46.4 10.5 116.5 

TABLE I: The position angles (8, <f>) and the orientation angles a of the four ground-based detectors (in units of degree). The 
angles (8, <j>) are given for a spherical coordinate on the Earth with the north pole at 8 — 0° . The orientation a is the angle 
between the local east direction and the bisecting line of two arms measured counter-clock wise. We use the abbreviations C, 
H, L and V as listed HHHl]. 

LIGO project and its advanced version (28|. Furthermore, the seismic noise level would generally become significant 
below ~ 10Hz for ground-based detectors. Therefore, we cannot simply apply the low-frequency-approximation for a 
world-wide-detector network with their typical noise curves. But it would be still interesting to examine the potential 
scientific advantage of fully using the four detector network in the long-term run. In addition our results might provide 
us with an useful insight for more realistic analysis. In this spirit, we numerically evaluated the generalized overlap 
functions C,n and Cvv- 

We obtained £77 = 0.07, a relatively low value compared with the results given in the previous subsection. For 
the circular polarization mode, our numerical results are £yy = —0.05 for CHL*V* , £yy = 0.011 for CH*LV* and 
£yy = —0.025 for CH*L*V with abbreviations C, H, L and V for the four detectors given in Table 1. The four 
detector network also has weak sensitivity to the circular polarization mode V2. 
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V. SIGNAL-TO-NOISE RATIO OF NON-GAUSSIANITY MEASUREMENT 

In this subsection we derive expressions for the signal-to-noise ratio of the non-Gaussianity measurement of a GW 
background made by superposition of burst signals. First, in subsection V.A, we write down the amplitude of the 
expectation values such as (C^iafc) and (K), including the fmiteness of the time segment T seg . Then, in subsection 
V.B, we estimate their fluctuations due to detector noises. 

A. signal strength 

We represent the time profile of a burst GW in the source frame by (X p (t), X c (t)), and take its Fourier transformation 
(k p (f),k c (f)) as follows 

/OO POO 
dtX p (t)exp[2TTift], k c (f)= / dtX e {t)exp[2irift], (53) 
-oc J — OO 

where an infinite time segment is used for above integrals. The expectation value for the total power of a burst is 
given by 

P(f) = (l\k P (f)\ 2 + Mf)\ 2 }s). (54) 

Our target background is a superposition of these bursts, and its spectrum is written by the burst rate r and the 
power P by (see e.g. (Tlj |) 

SgwU) = ^- (55) 

This spectrum is defined per solid angle and per polarization mode 0, Q . Here the well known normalized spectrum 
^Gw(f) is related to the spectrum Saw if) as 



■Vu !/i = 3 2 ^ /3 iW(/), (56) 



with the Hubble parameter Hq. In addition to the rate r and the power P, the burst duration Td is another principle 
parameter for characterizing the bursts and their background. 

Next we discuss analysis of the background using the Fourier modes (h p , h c ) transformed in a finite time segment 

Tseg as 

h p (f) = [ dtX p (t)exp[2mft], h c (f) = [ dtX c (t) exp[2mft]. (57) 

In these expressions, we omitted the label M for representing the initial time of the support of the integrals (see 
ea. (TlTJ|) V The expected number of bursts in the segment is given by 

q = rT seg max[l,T d /T seg ]. (58) 

Here the second factor in the right-hand-side is provided to include the effect that a single burst event is covered with 
multiple segments for Td > T seg . Meanwhile the amplitude I\ — (|/ip(/)| 2 + \h c (f)\ 2 ) in ea. (|4U)) is given by 

h - (\k P {f)\ 2 + |fc c (/)]| 2 )min[l,r se3 /T d ] = Pmin[l, T seg /T d ], (59) 

where the second factor represents the dilution of power due to a segment time T seg shorter than the signal duration 
T d . Using eas.(fI7 p (|iI J) ([58" J) and dHJ), we have 

(C 21ab ) = q (R a Rt) = = g?gg!gWfr (60) 

5 5 

Note that, with the spectrum Sqw defined in eq. (|55p . the amplitude for the traditional two point correlation (C^iab) 
does not depend on the burst duration Td- Nevertheless, this is not true for the Kurtosis parameter, as we see below. 
We can evaluate the expectation value (K) = q {RaRbR*R*d) in the same manner, and the result is given by 

(K) = ^rqlfW = ^ S GW PT seg Wwm[l, T seg /T d }. (61) 



12 




FIG. 3: Dependence of (K) and SNRk on the adjustable parameter T seg . The minimum value / _1 of the length T seg is set by 
the time resolution at the target frequency / in interest. We can estimate the signal duration Td by identifying the transition 
of the slopes for the signal (K) . 



Here the parameter W is defined in eq. (|48l) . and in the present case, we have 



Ir oih + (wV2 = oi (i\h P (f)\ 2 + ifac(/)i 2 ] 2 ) - cw ([h P (f)h c (fy - h c {f)h P urf) (r>| 

J l 2 ([l^(/)| 2 + l^c(/)| 2 ]) 2 

<n <DM/)I 2 + \kc(f)\ 2 ] 2 ) - Cw {[k P (f)k c (fr - k c {f)k p {fYY) 



((IM/)I 2 + IM/)I 2 )) 2 



(63) 



and it does not depend on T seg . Here we used ea.([59[) and the relations such as I2 — 
([|^p(/)| 2 + l^c(/)| 2 ] 2 ) {min[l, T se g/Td]} 2 ■ The segment time T seg is an adjustable parameter for signal analysis. 
When we increase the length T seg , there is a transition point at T seg = Td(> (A/) -1 ) where the signal (K) starts to 
decrease due to dilution of the power (see figure 3). We will discuss implication of this fact at the end of the next 
subsection. 

From eqs. (|55p and (|61[) we can derive a simple relation (K) oc t~ 1 S'q W for the event rate r and the background 
level Saw(f) (or Slow)- For a fixed amplitude Saw, the signal (K) becomes more Gaussian-like for a larger event 
rate r, as indicated by a smaller (K). This is quite reasonable from the central limit theorem. 



B. RMS fluctuations by detector noises 

In this subsection we evaluate the expected noise level for various statistical measures such as K or C-riab- As a 
brief summary of notations, we begin our study with providing the typical value of the optimal signal-to-noise ratio 
SNRsst for a single burst 

2(PA/)V 2 

SNRbSt = 5 1/2 5 V2 ' (64) 
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Here the factor 2 originates from the normalization associated with definition of signal-to-noise ratio [18| and the 
factor 1/5 1 / 2 is due to the averaging operation for the angular responses of detectors. Since we are interested in a 
GW background made by weak and undetectable bursts, we assume SNRsst < 1. 

The correlation C2iab is given by a summation of the signal products s^s^ for two detectors a and b over the 
Fourier modes i (see eq. (fT6|) ). Hereafter we assume that the GW spectrum Sew is smaller than the detector noise 
spectrum Sn (weak signal condition), and therefore the fluctuations for our statistical measures are dominated by 
detector noises. This corresponds to the condition when the correlation analysis becomes a powerful approach to 
detect a weak background buried among the detector noises, and is often assumed for theoretical analysis on GW 
backgrounds [f| (see also [f| for general cases). From eqs. ([55f and (|64|) we have Sgw/Sn = ^^SN R^ st r / /S.J and the 
condition for the assumption Saw "C Sn is given by 



r 



^NR L -«l. (65) 

The root-mean-square (rms) fluctuation for a product s a jS^ of each Fourier mode i is given by SNT seg / '(2\/2) (see 
eqs.(|10p and (fP2|) ) with the factor \j\pl for projection of the data toward the real axis on the complex plane [29] ]. 
Then the rms fluctuation of the second-order moment 5C 2 \ a b is given by 

5C 21ab ~^f, (66) 



where the factor 1/y/Nt (N t : number of Fourier modes) is the statistical suppression of fluctuation due to the 
summation of independent N t modes. 7 We can derive the same result for the fluctuation 8C22ab associated with the 
estimated moment C 22a b, and put 

5C 2 = SC 2 lab = SC 2 2ab- (67) 

For two aligned (jiab = 1) detectors, the expectation value of C 2 \ a b becomes (C 2 i a b) — s ' kSg ^ Ts " !! (see eq. (|6Tj)) ), 
and we obtain its signal-to-noise ratio as 

SNR C2 ^ <%^> = ^^(2T ofcs A / )V 2 . (68) 
6C 2 5 Sn 

The second power of this expression is essentially the same as the standard expression for correlation analysis in which 
we use the integral J df rather than the simple product x A/ [5j, [6J . 

From eq. (j2"0)) , the rms fluctuation of the parameter K due to the detector noises is estimated as follows 

SK - max[<SC 4 , C 2 SC 2 , {SC 2 ) 2 }. (69) 

Meanwhile the fluctuation for the term C4 is evaluated as 

C2 rp 2 

SC 4 ~ (70) 

The factors in this equation can be understood as in eq. (|6"6"|) . We can evaluate the ratios between the three elements 
in eq.([MD 8 

C 2 SC 2 16tt S GW < x (5C 2 f (71) 



Sd 5 S N Sd y/W t 

where we again assumed the weak signal case Saw "C Sn- If this condition does not hold, we need to deal with the 
contribution of the term C' 2 SC 2 . The relations (|7ip show that, in the present setting, the contribution from the term 
C4 dominates the fluctuation for the parameter K . Thus, for our weak signal case, we have 

«f~5C 4 ~Mk, (72) 



There can be a weak correlation between noises of different Fourier modes, especially with a short segment length T seg . But we neglect 
it here. 

The second expression shows that we can neglect the bias induced by the nonlinear combination for K in the weak signal condition. 
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and the signal-to-noise ratio for the parameter K is given by 

SNRk . W A ^f R r . (73, 
Afmax.[T d ,T seg \ 

Here we used the relation min[A, B] = 1/ max[A _1 , B^ 1 ]. In eq. ([42|) . the normalization factor 2/25 for the generalized 
overlap reduction functions (jj and (j/V is determined in order to simplify the pre-factor in eq. (|73[) . Since the time 
resolution at a frequency / effectively sets a limitation Td > / _1 ~ (A/) -1 , the minimum value of the denominator 
Afmax[T d ,T seg ] is — 1. 

From eq. ([75|) , we can understand that, even if the individual burst signal is too weak (SNRbst < 1) to be detected, 
its basic properties can be statistically studied with the amplification factor SNRc2- The factor SNRc2 is the 

_____ ^2 

signal-to-noise ratio for standard correlation analysis (eq. (|68[) ) and proportional to ^GwT obs ■ 

In relation to this, we compare magnitudes of SNR C 2 and SNRk- Since we consider a background made by weak 
undetectable bursts SNRb st < 1 and also have relations Sfmax[Td, T seg ) > 1 and W — 0(1), the signal-to-noise ratio 
SNRk would be generally smaller than that for the standard correlation analysis SNRc2 (see also [l]| for a better 
performance of a non-Gaussian statistic). 

As a function of the adjustable parameter T seg (> / _1 ), the signal-to-noise ratio SNRk is constant for a choice 
T S eg ^ Td, but it starts to decrease at T seg ~ Td (see figure 3). The expectation value (K) itself has a transition 
point at T seg ~ Td, as discussed after eq. (|61[) . Therefore, once we can detect the non-Gaussianity parameter K, we 
can estimate the typical burst duration Td by analyzing its background and identifying the transition. The result for 
SNRk also shows that it is statistically preferable to take a short segment time T seg , as assumed in this paper (see 
figure 3). 

In addition to the duration Td, we might also estimate the rate r and the characteristic power P of weak bursts, if 
(i) they are assumed to be the dominant source of the total GW background at a band and (ii) we can simultaneously 
detect signals C2 and K for the background. These two signals would provide us with two combinations proportional 
to rP and rP 2 . Then we can separately obtain the basic quantities P and r for the bursts by putting W ~ £77. 



C. GWs by Population III SNe 



In paper 1 we study a GW background caused by jet-like neutrino emissions at supernovae (SNe) of population III 
stars [22I f23j | . In the low frequency regime at / < 1Hz, the individual burst wave-form can be regarded as a simple 
step-function like time profile known as the burst of memory [301 ] . In relation to our analysis, its characteristic signal 
duration is Td ~ / _1 for a Fourier mode at a frequency /. The amplitude of the circular polarization mode V2 would 
be negligible, considering the linear emission pattern of neutrinos with I ^> |V| [30| . In this subsection, we reanalyze 
this background with the expressions presented so far. 

The followings are the characteristic model parameters extracted from [12] and used in Paper 1; the source redshift 
z ~ 15, the emitted neutrino energy E u ~ 10 55 erg, the mean anisotropy of the emission (q) ~ 0.03, and event rate 
r - O.Olsec" 1 . But we should notice that the actual values of these parameters are quite uncertain, since our current 
knowledge on the population III SNe is highly limited. For these model parameters, the amplitude of the background 
spectrum is given by View ~ 4 x 10~ 16 at / ~ 0.3Hz 11]. With the scaling relation flew °c rP, we characterize the 
burst background using the combination of parameters (Slew, 7 ") instead of the original ones (P,r). 

For the BBO noise spectrum in [2(|, the signal-to-noise ratio SNRbst oc {P-gw /r) 1 ^ 2 of the individual burst is given 
by 



SNR b , 




for / ~ 0.3Hz. Note that the relation ([55]) holds for the typical model parameters described above (see also figure 
4). We have the signal-to- noise ratio SNR c2 for the correlation analysis with two overlapped detectors (i.e. 7/ = 1 
as the A1-A2 pair in subsection IV) 

From eq. ([75)) the signal-to- noise ratio SNRk for the Kurtosis parameter K becomes 



SNRk - W SNR 2 bst SNR C2 , 



(76) 
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FIG. 4: Analysis for gravitational wave background from population III SNe with BBO. We have the scaling relations SNRc2 oc 
rSNRl at and SNRk oc rSNRf at . The typical model parameters used in the main text are SNR^st = 0.6 (individual burst 
strength) and r = O.Olsec -1 (burst rate). We put the observational time T f, s = lOyr, the bandwidth A/=0.3Hz and W = 5/14. 
The dashed line at the upper right is Sgw = Sn for the weak signal condition used in our formulation (see ea. (!65|l ). The 
vertical dashed line SNRf, st = 1 shows the boundary for weak undetectable signals. 



for the optimal setting T seg ~ / _1 of the segment time. 

With the generalized overlap reduction function Qj — 5/14 (see eq. ([50|) ) for BBO and the relation W ~ O1I2/I1 ~ 
C/i 9 for the non-dimensional parameter W defined in eq.(|55)l. we obtain the following result that is identical to eq.(18) 
in Paper 1 

SNRk ~ ±SNHl, SNHc, ~ 10 (j^)" (^)" /2 (jjjj^r)-' ■ (77) 

This result shows that, while the relevant astronomical parameters are highly uncertain, the non-Gaussianity signature 
(K) might be detected for the population III SN background with BBO. However, if the combination of the background 
parameters fi^^r -1 is smaller than ~ 2 x 10 _29 sec, the detection would be difficult. In figure 4, we show the schematic 
picture of the signal-to-noise ratios SNRc2 and SNRk as functions of SNR^st and r. 



VI. SUMMARY AND DISCUSSIONS 



In this paper we discussed an observational method to analyze a GW background made by a superposition of weak 
burst signals that are individually undetectable with small amplitudes. Our approach is to measure non-Gaussianity 
of the background induced by the discreteness of the bursts. The basic strategy was briefly introduced in Paper 1. 
This paper is a follow-up study with a significant enhancement on the basic formulation to clarify how we can study 
polarization properties of burst signals in response to the geometry of a detector network. 

We find that our method enables us to measure linear combinations of two interesting parameters I2 and V2 defined 
by averaged squares of the Stokes parameters I and V of individual burst signals. These two parameters I2 and V2 
can be separately estimated by devising data analysis with using multiple detectors, and they might provide us with 



Note that the ratio I2/I1 depends on the probability distribution of the burst amplitudes. The numerator I2 is more affected by the 
stronger (but undetectable) ones. 
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an important insight to discriminate a GW background once detected by the standard correlation analysis with two 
detectors. 

In this paper, the shortness of the signal duration has been very important to show up the discrete nature of 
bursts by localizing their power in a limited time domain. As an extension of our method, it would be interesting to 
study non-Gaussianity induced by almost monochromatic GW sources. In this case the individual signal is localized in 
the frequency domain, and it would be advantageous to use a relatively long time segment for Fourier transformation. 

We have made several assumptions and simplifications on data analysis. In order to apply the present method to 
real data of detectors, further careful works should be done. These include potential effects of non-Gaussianity and 
correlation of detector noises, and implementation on the short-term Fourier transformation. For example, if detectors 
have correlated non-Gaussian noises, performance of our method would be limited. In particular, two data streams 
(Ai,Ei) obtained from one set of BBO (or DECIGO) could have correlated noises, depending on the symmetry of 
the three vertexes of the unit. This might reversely provide interesting requirements on the mission designs of these 
future projects. 

The author would like to thank T. Tanaka for helpful discussions. He also thanks anonymous referees for invaluable 
comments to improve the draft. This work was supported by Grants-in-Aid for Scientific Research of the Japanese 
Ministry of Education, Culture, Sports, Science, and Technology 20740151. 



APPENDIX A: FOURIER TRANSFORMATION WITH A FINITE TIME SEGMENT 

Here we summarize basic correlation properties of Fourier modes transformed with a finite time segment T seg . We 
first express a time-domain signal a(t) using the continuous Fourier modes A(f) (defined with infinite time segment) 
as follows 

/oo 
e- 2 ^A(f)df, (Al) 

and 

A(f) = / e^*o(t)4f. (A2) 



We assume that the signal a(t) is a real function, or equivalently A(—f) = A(f)*. We define the power spectrum 
S(f) of the modes A(f) by 

(A(f)A(fT) = ~S(f)5(f-n (A3) 
Next we evaluate the Fourier modes A{f,T seg ) defined with a finite segment length T seg 

A(f,T seg )= [ Tse9 eWaftdt (A4) 
Jo 

at discretized frequencies / = N t x T"* (N t : integer). From eq. (|Al|) we can easily show the following relation; 

A(f,T seg ) = f S ' 3 e^a{t)dt = f°° df ^ W ~ A(f')e^-^ T «' ■ (A5) 

Taking into account the finiteness of the segment T seg (see eq. (fT2|) ). the power spectrum S(f,T seg ) for the present 
case is given by 

S(f,T seg ) = ±- (A(f,T seg )A(f,T seg )*) = r d f' s ^Ml-I^As(n. (A6) 

J -co [7T(/ - f')\ T seg 

With the following asymptotic relation for the Dirac's delta function 5(x) 

t^H^= s[i) - (A7) 



17 



the estimated power spectrum £>(/, T seg ) coincides with the original one S(f) in the limit T seg — > oo as 

lim S(f,T seg ) = S(f). (A8) 

Eq. (|A6[) shows that, for a finite T seg , the spectrum S(f, T seg ) has a dominant contribution around the target frequency 
/ with a bandwidth A/ ~ T~ g . 

Now we evaluate the correlation 2 (A(f, T seg )A(f, T seg )) /T seg corresponding to the moment G2 in eq.([J|. From 
eq. (|A5|) , we have 

J- { A(f,T seg )A(f,T seg)) = * r r HAf ~ f J;^ [J + n (A9) 

1 seg 1 seg J — oo " \J J ) 

In the limit T seg — > 00 we have 

2 IP 00 
lim — (A{f, T seg )A(f, T seg )) = lim — / d/'[<5(/ - /') + J(/ + /')]S(/')e 2 " /T - = 0, (A10) 

Jseg > oo J- seg 1 seg >oo 1 seg J — oo 

and the expectation value for the amplitudes of real and imaginary parts of the mode A(f,T seg ) are the same, as 
shown from the real part of ea. l|A10p . In addition, the imaginary part of ea. (|A10p represents that the correlation of 
the two parts vanishes. However, these properties do not hold at a finite T seg - Therefore, in contrast to the limit 
Tseg — > oo, we need to keep the correlation (A(f,T seg )A(f,T seg )) for a finite T seg - 



APPENDIX B: ANALYSIS IN PAPER 1 



In this appendix, we comment on the two overlooked points in Paper 1. They are caused by improper handling of 
Fourier modes derived with a finite time segment T seg , as discussed in Appendix A. 

In Paper 1, we consider the simple product (s_b 1 s^ 2 ) 2 for analyzing non-Gaussianity of a burst GW background. 
Here the data sex — Hex + nEX and SA2 = Ha2 + n^i are made by GW signals Hex and Hai and detector noises 
Hex and nA2 with the labels Ai and E\ for detectors defined in subsection IV. A. But, we have the correlations for 
GW signals (HE 1 (f)HE 1 {f)) ^ and (H A 2 (f)* 'H A 2 (/)*) whose amplitudes can become large for a short segment 
T S eg- We thus need to subtract the relevant terms in order to extract a non-Gaussianity signature, as outlined in 
section II. A. 

Furthermore, with a finite T seg we also have the correlations of detector noises (nE 1 {f)*nE 1 (f)*) ^ and 
(nA 2 {f)* n A 2 (f)*) f° r tne simple product {se 1 s* A2 ) 2 , and the detector noises have a contribution to the ex- 
pectation value of the product (sb 1 s^ 2 ) 2 . This qualitatively changes the statistical character of the present problem 
and decreases the signal-to-noise ratio of the non-Gaussianity measurement, compared with the combination, such as 
( s e 1 se 2 s *a 1 s a 2 ) use d in subsection IV. A. 
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